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ABSTRACT 


The asymptotic solutions, described by Tennekes and Lumley (1972), for surface flows in a channel, 
pipe or boundary layer at large Reynolds numbers are revisited. These solutions can be extended 
to more complex flows such as the flows with various pressure gradients, zero wall stress and 
rough surfaces, etc. In computational fluid dynamics (CFD), these solutions can be used as the 
boundary conditions to bridge the near-wall region of turbulent flows so that there is no need to 
have the fine grids near the wall unless the near-wall flow structures are required to resolve. These 
solutions are referred to as the wall functions. Furthermore, a generalized and unified law of the 
wall which is valid for whole surface layer (including viscous sublayer, buffer layer and inertial 
sublayer) is analytically constructed. The generalized law of the wall shows that the effect of both 
adverse and favorable pressure gradients on the surface flow is very significant. Such an unified 
wall function will be useful not only in deriving analytic expressions for surface flow properties but 
also bringing a great convenience for CFD methods to place accurate boundary conditions at any 
location away from the wall. The extended wall functions introduced in this paper can be used 
for complex flows with acceleration, deceleration, separation, recirculation and rough surfaces. 

1 INTRODUCTION 


An asymptotic solution for the inertial sublayer in a channel or pipe flow at large Reynolds 
numbers can be written as (Millikan, 1938) 



where JJ is the mean velocity, u T is the skin friction velocity defined by the wall stress r w as 
u T = y/\r w \/p, y is the normal distance from the wall, v and p are the viscosity and density of 
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the fluid, k « 0.41 and C ~ 5.0. Eq.(l) is also theoretically valid and only valid for a flat plate 
boundary layer, but it has been applied to other wall bounded flows with some successes despite 
its formal validity. For a boundary layer with an adverse pressure gradient and zero wall stress, 
Tennekes and Lumley (1972) derived another asymptotic solution which reads 



where u p is defined by the adverse wall pressure gradient as u p = [(i//p)\dP w /dx\} 1 ^ 3 , and a « 5, 
p « 8 according to the experimental data of Stratford (1959). Eq.(2) has not been paid much 
attention in computational fluid dynamics. Apparently, Eq.(l) will become erroneous for flows 
near separation or re-attachment points because there the wall stress, hence the skin friction 
velocity, is nearly zero. On the other hand, Eq.(2) will not be valid for boundary layer flows with 
a small or zero pressure gradient because u p is nearly zero. 

In this paper, we will briefly repeat the analyses of Tennekes and Lumley and introduce a more 
general asymptotic solution for the surface flow valid for both the zero or nonzero wall stress and 
the zero or nonzero wall pressure gradient. Therefore, the solution can be used for flows with 
acceleration, deceleration, separation and recirculation. 

The basic idea is to assume, at large Reynolds numbers, the existence of a surface layer distinct 
from the outer layer in a boundary layer flow. The existence of the law of the wall in the surface 
layer and the existence of the velocity-defect law in the outer layer will lead to an asymptotic 
solution for the surface flow in the region where 

| < 1 , f - » 1. (3) 

where, 6 represents the thickness of the boundary layer, t v is the length scale related to the 
viscosity of the fluid which will be defined later. The region in which Eq.(3) holds is called the 
inertial sublayer. In the vicinity of the wall where y/t v is of order one, the turbulent stress is 
significantly suppressed and this region is called viscous sublayer. An asymptotic solution can 
be obtain for both the inertial sublayer and the viscous sublayer. The region between these 
two sublayers is called buffer layer where the turbulent and viscous stresses are of same order. 
A simple model of turbulent stress leads to an expression which matches both the viscous and 
inertial sublayers and leads to an unified law of the wall, similar to the one proposed by Spalding 
(1961). We will start first with the boundary layer flow over a smooth wall. A general asymptotic 
solution for the surface layer will be obtained. The effect of rough surfaces on the solution will 
then be considered. 


2 WALL BOUNDED TURBULENT FLOWS 


The equations of motion for steady two-dimensional incompressible flow in a Cartesian coordinate 
system are 


dU_ dV_ 
dx + dy 


0 , 


( 4 ) 
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( 5 ) 


dU dU _ 1 dP duu duv d 2 U d 2 U 

dx dy p dx dx dy dx 2 dy 2 ’ 

dV 1 dP duv dvv d 2 V d 2 V 

dx + dy p dy dx dy dx 2 dy 2 

For simplicity, we are considering boundary layer flows in which the effect of wall curvature 
is neglected. The existence of the wall and the no- slip condition at the wall will create, at a 
sufficiently large Reynolds number, a very thin surface layer near the wall, which is distinct from 
the outer layer of a turbulent boundary layer. In the surface layer, the flow is largely affected by 
the viscosity, the governing equations (4-6) can be significantly simplified and their solutions are 
of certain forms called the law of the wall. On the other hand, in the outer layer of a boundary 
layer remote from the surface layer, the flow is less or nearly not affected by the viscosity, the 
solution of the above equations is of another forms called the velocity- defect law. 


2.1 The law of the wall 


In a coordinate system with the wall at y = 0, the boundary layer flow in the half-plane y > 0 
is mainly in the x direction, except near the separation or re-attachment point. In general, the 
wall stress t w and the wall pressure gradient ( dP/dx) w are non-zero and can be either positive 
or negative with respect to x direction. We may define a velocity scale u c using these two wall 
parameters as follows, 



Thus defined u c will never become zero in any boundary layer flows with either zero wall stress 
or zero pressure gradient because u T and u p cannot be zero at the same time. With u c we may 
define a viscous length scale t v = v/u c . This viscous length scale is usually very small comparing 
with other length scales of the boundary layer, for example, the boundary layer thickness 6 and 
the downstream length scale L. That is, u c 6/v >> 1 and u c Lj v 1. 


Let us start with the simplification of the governing equations (4) -(6) for the flows in the thin 
surface layer using an order of magnitude analysis (see Tennekes and Lumley, 1972). We use u c 
to scale both the mean velocity U and turbulent velocities u, v. Let L be the downstream length 
scale and i v — v/u c be the length scale in y direction. With dU/dx ~ uJL and dV/dy ~ Vji v , 
the continuity equation (4) gives V ~ u c l u f L. The left hand side of Eq.(6) is then of order 
u\lvlL 2 . The orders of magnitude of the turbulent terms in Eq.(6) are 


dvv 

dy 


= 0 


the viscous terms in Eq.(6) are of order 

d 2 V 



duv_ n (K' 

dx \L , 


v 


d 2 v 

dx 2 


Ut il 


o 

\L L 2 , 


( 8 ) 

( 9 ) 
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Because u c Ljv > 1 or < 1, the major turbulence term, dvv/dy , must be balanced by the 
pressure term in the surface layer, that is 


1 dP_ dm 

p dy + dy 


( 10 ) 


Integration of Eq.(10) from the wall to y (within the surface layer, i.e., y/6 < 1) and differentiation 
with respect to x lead to an expression for dP/dx: 


1 dP _ 1 dP w _ dm 
p dx p dx dx 


( 11 ) 


where, P w is the pressure at the wall. With Eq.(ll), the pressure term in the x-momentum 
equation (5) can be expressed in terms of the wall pressure. 


Now we may estimate the various terms in Eq.(5) in the surface layer as follows 



The major terms in Eq.(12) are duv/dy and vd 2 U/dy 2 if l v jL <C 1. Therefore, the x-momentum 
equation (5) can be approximated in the surface layer as follows 


8 ( _ dU\ 1 dP w 

Tv\- mJrV W = -p-^- 


(13) 


Integration of this equation from y = 0 to y (within the surface layer) yields 

. dU t w y dP w 

- UV + v— = — -(- — — , 
dy p p dx 


(14) 


where r w = p(dU/dy) y= o is the wall stress. The no-slip condition at the wall for turbulent 
velocities has been imposed. This is a general equation for surface flows under the condition 
u c Ljv > 1. Eq.(14) indicates that, in general, the surface flows are affected by both the wall 
stress and the wall pressure gradient. The relative importance of these two terms depends on the 
flow situation. For example, for a boundary layer flow with a strong adverse pressure gradient, 
the wall stress t w can become very small and even vanish. In this case, the adverse wall pressure 
gradient controls the surface flow and the total shear stress (the left hand side of Eq.(14)) is not 
constant across the surface layer. On the other hand, if the pressure gradient is zero or small 
compared to the wall stress, then the wall stress dominates the surface flow and the total stress 
is constant or nearly constant across the whole surface layer. 


The flows we want to consider here include acceleration, deceleration and even recirculation. 
Therefore, the two terms on the right hand side of Eq.(14) could be of the same order of magnitude, 
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or one is large than the other. Note that because Eq.(14) is linear, we may deal with these two 
factors separately by decompose U and —uv into two parts (Tennekes, 1968, has shown with a 
multivariate asymptotic technique that this is a valid procedure). Following Tennekes and Lumley, 
we write 


U = Ui + u 2 . 

(15) 

— UV = —(uv) 1 — (uv) 2 - 

(16) 


The first part, represented by U\ and — (m?)i, is associated with the wall stress r w / p only and the 
second part, represented by U 2 and —(uv) 2 , is solely related to the pressure gradient (y/p)dP w /dx, 
i.e., 


v , dUi T w 
-(«»), + dy - p , 

(17) 

( A , dU 2 y dP w 

(uv) 2 + v = . 

dy p ax 

(18) 

The first part of the flow in Eq.(17) has only one characteristic velocity 

u T and one characteristic 

length v/u T . Similarly, the second part of flow in Eq.(18) has only one 

characteristic velocity u p 

and one characteristic length u/u v . Therefor, the nondimensional form of Eq,(17) and Eq.(18) 

can be written as 


~(uv)l d(Ui/u r ) _ T w 

u\ d(u T y[v) pul' 

(19) 

~(uv) 2 d(U 2 /u p ) _ ( u p y\ v dP w j dx 
u 2 p d(u v y/u) ~ V v ) p 

(20) 

There are no additional parameters in the boundary conditions on Eq.(19) and Eq.(20) because 
the boundary conditions are homogeneous (both the mean velocity and the turbulent stress are 
zero at y = 0). Therefor, the solution of Eq.(19) and Eq.(20) must be of the forms: 

a_^ A (3a). 

u T pu * \ v J 

(21) 

— (tu;)i _ t w / u T y\ 

u 2 ~ pu 2 . 91 V v ) ' 

(22) 

and 


U 2 i/dP w /dXj fu p y\ 
U p p U 3 2 V V J ’ 

(23) 

(uv) 2 is dP w /dx / u p y\ 

^ ~ P Up 92 V v J ' 

(24) 


Eq.(21)-Eq.(24) are called the law of the wall. 
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2.2 The velocity-defect law 


In the outer layer of a boundary layer, the boundary layer thickness 8 is the only appropriate 
length scale. If u r is the only characteristic velocity for the first part of the flow, then the first 
part of the velocity-defect (U\ - Uq)I{t w / pv? r ) is a function of y, 8, u T only, where Uo is the velocity 
at the edge of the boundary. Therefore, we have a system of four quantities with two dimensions. 
From II theory of dimensional analysis, only two independent normalized quantities can be formed 
and we may write this system as 


or 


f UUi-Uo)/»t 

l r w jpu\ 


Ul -U 0 T w 

= 2' tl 

u T pu‘ 



(25) 


For the second part of flow U 2 , we may use the above same argument to obtain a relation for U 2 


U 2 ~ Uo _ v_ dPw[dx_ f y\ 

u p p uj 2 \8/ ' 


The relations (25) and (26) are called the velocity-defect law. 


3 ASYMPTOTIC SOLUTIONS FOR SURFACE LAYER 


3.1 The inertial sublayer 


If the Reynolds number, Re = u c 8jv , is large enough, an overlapping layer between the surface 
layer and the outer layer may be developed. This overlapping layer is called inertial sublayer, 
where y/8 1 and u c y/i/ > 1. The existence of the inertial sublayer can be easily seen from the 

following relation 

(27) 


V 

8 


u c y v 


v u c 


For example, if u c 8jv > 10 4 , there will exist a region where u c y/v > 10 2 and y/8 < 10 -2 . In the 
inertial sublayer, the law of the wall in the surface layer should match the velocity- defect law in 
the outer layer. Following Tennekes and Lumley to equate the mean velocity gradient dUi/dy 
calculated from Eq.(21) and Eq.(25), we obtain 


f u T y \ df 1 f y\ dR\ 

\ v ) d(u T y[v ) d(y/8)' 


(28) 


Eq.(28) indicates that the both sides must be equal to a same constant, 1 /k, say. And this leads 
to a logarithmic velocity profile in the inertial sublayer: 


Ui 

u r 


pv? T L« \ v J 


+ C . 


(29) 
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To match the case with zero pressure gradient boundary layer, k = 0.41, C = 5.0. Eq.(29) can be 
rearranged using u c as 

= (30) 

U c pu~ L KU C \ V J 

where 


= 


U r 


U c 


-ln(-) 

K \U C J 


+ C 


( 31 ) 


C\ is shown in Figure 1, it goes to zero as u T vanishes (u p /u c — >• 1.0). 


Using the above same argument for the second part of the flow U 2 , based on the law of the wall 
(23) and the velocity- defect law (26), we may obtain the solution: 


U 2 _ v dP w / dx 

Up p u ^ 


K“M- 


(32) 


To match the case with zero wall stress boundary layer, a = 5.0 and /? = 8.0. Eq.(32) can be 
rearranged as 


U 2 ^ dPm j dx 


U r 


LWm') 

. U c \ v J 


+ C2 


where 


C 2 = — fain ^ + P 

u c L \u c J J 

C 2 is shown in Figure 1, it goes to zero as u p vanishes. 


(33) 


(34) 



Figure 1: Coefficients C\ and C 2 in the law of the wall 
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( 35 ) 


Finally, the total velocity U = Ui + U 2 can be written as 


u _ 

7~vj 

liillnfMl 

+ C\ 

vdP w jdx 

"I" 0 

[ q Wm) + c J 

Uc 

■ pul 

.* U c \ V ) 


p u f 

[ u c \ v J J 


With Eq.(35) and Eq.(14), it can be shown that, in the inertial sublayer ( u c yjv >1), the total 


turbulent stress is 


— uv = 


Tw y dPyj 

p pdx 


(36) 


uv t w f u T \ dP w /dx f Up \ UcV f37) 

u l pu % \u c J \dP w /dx\ \n c / v 

Equations (35) and (36) are the asymptotic solutions for the surface flow in the inertial sublayer, 
where u c y/u > 1 and y/6 < 1. These equations are called wall functions in CFD, because these 
relations can be used as the near-wall boundary conditions for CFD calculations. The effect of 
pressure gradients on the mean velocity profile described by Eq. (35) is shown in Figure 2. The 
left of Figure 2 shows the strong effect of adverse pressure gradients on the mean velocity profile 
as the ratio of u p /u c varies from zero to one (which corresponds respectively to the zero pressure 
gradient boundary layer and the boundary layer about to separate). In the boundary layers with 
favorable pressure gradients, the skin friction velocity u T increases with the pressure gradient so 
that u p /u c = u p /(u T + Up) is always less than one and it is usually a small value. For example, in 

a fully developed channel flow, it can be shown that u p /u c = 1/(1 + R l Jh), where the Reynolds 
number is defined as R e h = U ma xh/v, h and U max are the half width of the channel and the center 
line mean velocity, respectively. Therefore, if R e k > 10 4 , then u p /u c < 0.044, hence, the effect of 
pressure gradients on channel or pipe flows will not be very significant if the Reynolds number 
R e h is sufficiently large. However, for accelerating boundary layer flows, the value of u p [ u c could 
be much larger than 0.05 before the turbulence has been suppressed by the acceleration. The 
right figure of Figure 2 shows the effect of favorable pressure gradients on the mean flow up to 
u p /u c = 0.3, which corresponds to the flows with extremely large favorable pressure gradients. 
The effect is significant. The effect of pressure gradients on the turbulent stress, described by Eq. 



Figure 2: Effect of pressure gradients on the mean velocity 
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(36) or Eq. (37), is shown in Figure 3. Again, both the adverse and favorable pressure gradients 
have the strong effect on the turbulent shear stress and there is no constant shear stress layer for 
boundary layer flows with large pressure gradients. The constant shear stress layer only exists for 
flows with zero pressure gradients, or u p /u c << 1. In addition, we see —Wv changes its sign for 
some large favorable pressure gradients in the right of Figure 3, which may indicate the limitation 
of Eq. (36) or Eq. (37). 



CM O 

^ 0.5 

3 
V 


Favorable pressure gradient 


— 


\ 


" ■ u^/u. e 0.0 

— - u^/u, * 0.05 

— u^u. B 0.1 

— — u^/u. * 0.2 

\ i ij\x t ■ 0.3 

\ 


\ 


uj 


u c y/v 


200 300 400 500 


Figure 3: Effect of pressure gradients on the turbulent stress 


3.2 The viscous sublayer 


In the vicinity of the wall, for example, u c y/v < 5, the viscous effect dominates the flow. The 
turbulent stress -W is vanishingly small compared to the viscous stress udU/dy even though the 
flow is still quite disturbed and is not laminar. In fact, the ratio of the turbulent rms velocity 
u' to the mean velocity U at the wall is finite, i.e., (u'/U) y -+ o = const.. In this so-called viscous 
sublayer, the turbulent stress can be neglected and the velocity profile can be obtained from 
Eq.(14) as follows 

TT /. .A 1 J 7"> /J_ /-. *.\ 4 

(38) 


U_ _ Tw_ 1 v dP w /dx / u c y Y 

u c pu 2 \ v ) 2 p uf \ v j 


or, 


Twy ivdPw fy \ 2 

p v 2 p dx \vj 


(39) 


The effect of pressure gradients on the flow in the viscous sublayer is shown in Figure 4. Again, 
the effect of pressure gradients is significant. 
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Figure 4: Effect of adverse pressure gradients on the flow in the viscous sublayer. 

4 EFFECT OF ROUGH SURFACES 


Now let us consider the effect of rough surfaces on the surface flow. Let h denote an rms roughness 
height. If the ratio h/6 is small, then the roughness will not affect the velocity-defect law. However, 
the law of the wall in the surface layer may need modification. Let us define y = 0 as the average 
vertical position over the rough surface. Apparently, we are interested in the flow only at y > h, 
because at y = 0 the flow field is not even defined. Now, the surface layer over a rough surface 
has another length scale h, in addition to v/u c . The ratio of the two is the roughness Reynolds 
number Rh = u c hjv. If Rh is of order one, Rh < 5 say, one may expect that the roughness will 
have no effect on the surface flow because the roughness elements are submerged in the viscous 
sublayer where no turbulent stress can be generated, even though much the flow is disturbed. It 
can be shown that Eq.(14) is still valid for the surface flow with rough surfaces as long as h/6 < 1 
and u c L/v >• 1, h/L <C 1. Following Tennekes and Lumley’s argument^, the law of the wall can 
be written as 


Ei = jw A to , 

U T pUj \ v J 

(40) 

U 2 vdP w /dx / u p y ^ \ 
u p P Up 2 V v ) ’ 

(41) 

u T pu l T \h ) 

(42) 

v 1 = t dpjd* j, y 

Up p u* \h J 

(43) 


The velocity- defect law in the outer layer of a boundary, Eq.(25) and Eq.(26), will be independent 
of the surface roughness as long as h/6 *C 1. The matching of the velocity derivative in the inertial 
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sublayer will lead to the logarithmic velocity profile with an additive function depending on Rh- 


U_ 

u c 


L 




K U c 


u ) 


+ Ci(Rh) 


v dP w /dx 




a^lnf^)+C 2 (^)l, (44) 

u c \ v J J 


or 


£ 

U c 


pul 


-—to (t)+C'i(Rh) 


K U c 


v dP w /dx 


U'i 


°vA£) 


+ C' 2 {R h ) 


(45) 


For Rh < 5, as mentioned earlier, the surface roughness is too small to affect the surface flow. 
Hence, in Eq.(44), C\{Rh ) C\ , C2(Rh) C 2 when Rh < 5. For Rh > 5, Ci(Rh) an d C 2 (Rh) 
may depend on the roughness Reynolds number Rh* The detailed relations must be determined 
by experiments. 

However, for large values of Rh (> 30, say), Tennekes and Lumley^ have shown that the co- 
efficients C[{R h ) and C' 2 (R h ) in Eq.(45) should be independent of R h . If these coefficients are 
absorbed in the definition of h, then Eq.(45) may be written as 


U_ 

u c 



v dP w j dx 



(46) 


With Eq.(44) or Eq.(46) and Eq.(14), it can be shown that the turbulent stress in the inertial 
sublayer (u c y/ v > 1) is 

y dP w 


- uv = — H 

p p dx 


(47) 


Now let us consider whether there exists a viscous sublayer on a rough surface. Apparently, for the 
case of large Rh, the roughness elements are well above the viscous sublayer defined by u c y/v < 5. 
These elements generate turbulent wakes and are responsible for essentially inviscid drag on the 
surface. Therefore, no viscous sublayer exists and Eq.(46) may be used down to the “wall”, where 
y = h,U = 0. 

On the other hand, however, for small Rh < 5, the viscous sublayer may exist. It can be shown 
that Eq.(39) is valid with the following modification: 

U = ^{y - h) + \^^-{y - hf. (48) 

pv K 2 pv 

where U = 0 at y = h. Eq.(46) and Eq.(48) indicate that the “effective wall” for rough surfaces 
is at y = h. 


5 Unified wall function 


We have discussed the asymptotic solutions for surface flows in the inertia sublayer and the viscous 
sublayer. The region between these two sublayers, i.e., 5 < u c yjv < 30, is called buffer layer. 
In this region, the viscous and turbulent stresses are of same order. No theoretical asymptotic 
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solution can be obtained for this layer. However, by using a simple model for the turbulent stress 
in the buffer layer and a proper matching procedure, we are able to obtain a single analytic 
function, unified wall function, for the whole surface layer which includes viscous sublayer, buffer 
layer and inertia sublayer. Such an unified wall function (without considering the effect of pressure 
gradients) was first constructed by Spalding (1961). Here, we will construct an unified wall 
function with the effect of pressure gradients so that it can be used for flows with acceleration, 
deceleration, separation and recirculation. 


5.1 Unified law of the wall 


Let us start with Eqs. (17) and (18). In the buffer layer, -(uv)i and v ^ are assumed of same 
order. The same is true for — (i uv ) 2 and v^-. The buffer layer is considered very close to the 
wall. Therefore, we may assume that the following near-wall behaviors are valid throughout the 
buffer layer, 

— (uu)i ~ y 3 + 0(y 4 ) , U\ ~ y + 0(y 2 ) , hence - ( uv)i ~ Uf + 0 (llfj (49) 

~(uv) 2 ~ y 3 + 0(y 4 ) , U 2 ~ y 2 + 0(y 3 ) , hence - (uv) 2 ~ U^ 2 + 0 (u 2 /2 ) (50) 


If we further model -(uu)i and ~(uv) 2 with an eddy viscosity concept that they are proportional 
to the mean velocity gradient, and note that dU\/dy is order of one and dU 2 /dy is of order y or 
U^ 2 , then we may write 


- (uu)i oc v 




dU i 
dy ’ 


— (uv) 2 OC V 




dU2 

dy 


(51) 


With theses models, Eqs. (17) and (18), in the region from the wall to the buffer layer, can be 
written as 


v 


1+C'S + 0 |-L 


1 + C"— + o 

ViT) 



dlh = —dy 
P 

jtt y j 

du 2 = -^— d y 

p ax 


(52) 

(53) 


Integrate the above equations, we obtain 


where 


y+ = u+ + c'u + 4 + o (f/+ 5 ) , 

(y +) 2 = u+ + c"uf + o {yf) 

y+ _ U rV y+ _ M Ul u+ _ 2 j dpw/dx\ U 2 

T v ’ p v ’ 1 t w u T 2 dp w /dx Up 


(54) 

(55) 

(56) 


Equations (54) and (55) describe the behavior of U\ and U 2 in the buffer layer which matches 
their behaviors in the viscous sublayer. The coefficients in the higher order terms, C', C", • • • , are 
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unknown. However, they will automatically be determined by matching Equations (54) and (55) 
with the asymptotic solutions in the inertial sublayer. In the inertia sublayer, the asymptotic 
solutions (29) and (32) can be written as 

Yj = exp(-nC) exp(nU^) (57) 

(y+) 2 = exp(-2a/0 ) expiUf /a) (58) 

Their series expansions in terms of U+ and are 

Y+ = exp(-KC) 1 + nUf + ^ ( kU f ) + ^ ) +‘” 

( y +) 2 = exp(-20/a) 1 + Z#/a + • • •] 

Now, in order to have an expression for U+ which behaves like Eq. (57) when u T y/v is large (say, 
> 30), otherwise, Eq. (54) when u T yjv < 30, the following form is the simplest: 

Y+ = Ui + exp(-nC) exp(nU + ) - 1 - kU+ - - ) “ g ( K ^i + ) ( 5 ®) 

This analytical expression has been supported by many experimental data and is referred to as 
the Spalding’s law of the wall. Similarly, by matching Eqs. (58) and (55) we obtain another 
analytical expression, the law of the wall for flows with zero wall-stress: 

(Y p + ) 2 = Ut + exp(—20/a) [exp(U+/a) - 1 - U+/a] (60) 

The unified laws of the wall for zero pressure gradient and zero wall stress, (59) and (60), are 
shown in Figure 5. 



Figure 5: Unified laws of the wall 


From the application point of view, we need the inverse form of Eqs. (59) and (60). Unfortunately, 
it is impossible to obtain a single analytical inverse form for both (59) and (60). Therefore, we 
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suggest the following piecemeal functions: 


Ui 

U T 

£2 

U p 




dp w /dx 

\dp w /dx\ 


MY*) 


where f\ and / 2 are the piecemeal fitting functions defined as follows 


fl(Y r + ) = { 


f2(Y p + ) = \ 


a 1 y+ + a 2 (y +) 2 + a 3 (y T + ) 3 

60 + hY T + + b 2 (y+) 2 + b 3 (y+) 3 + 64 (y r + ) 4 

c 0 + Cl y+ + c 2 (y+) 2 + c 3 (y+) 3 + c 4 (y+) 4 

iio g (y+) + c 

K 

a 2 (y p + ) 2 + a 3 (y/) 3 2 

6 0 + ^y p + + *2 (y/) 2 + h (y+) 3 + b 4 (y+) 4 

c 0 + Cl y+ + c 2 (y+) 2 + c 3 (y+) 3 + c 4 (y+) 4 
a log (y+) +8 


if Y+ < 5; 
if 5 < Y+ < 30; 
if 30 < y+ < 140 

if 140 < y+ 

ify+ <4; 
if 4 < y p + < 15; 
if 15 < y+ < 30; 
if 30 < y+ 


The coefficients in /1 


fll 





1.0 

1.03-02 

-2.93 - 03 



bo 

b 1 

b 2 

&3 

b 4 

-0.872 

1.465 

-7.02 3 - 02 

1.663 - 03 

-1.4953 - 05 

Co 

Cl 

C2 

C3 

c 4 

8.6 

0.1864 

-2.0063 - 03 

1.1443 - 05 

-2.5513 — 08 


The coefficients in / 2 


0*2 

&3 




0.5 

-7.313-03 




bo 

hi 

h 2 

b 3 

b 4 

-15.138 

8.4688 

-0.81976 

3.72923 - 02 

-6.38663-04 

co 

Cl 

C2 

C3 

c 4 

11.925 

0.93400 

-2.78053 - 02 

4.62623-04 

-3.14423 - 06 


The inverse forms of Eq. (63) and Eq. (64), are also shown in Figure 5. 





5.2 Unified wall function 


Now, we may construct the unified wall function for a general surface flow using the results 
described in the previous section. Since 

U _ Ui + U2 _ Ut Ui + UpUi 

Hq vl^ u c y>T u c Up 


we may write 


u_ 

Uc 


T w U T £ 

2 * ■ 1 

pU% U c 



dp w /dx Up f + Up\ 
\dp w /dx\ u c J2 \ «c/ 


(65) 


where 



The left and right figures in Figure 6 show respectively the effect of adverse and favorable pressure 
gradients on the surface flow. As the adverse pressure gradient increases, the ratio of u p /u c 
increase from zero to one and u T /u c decreases from one to zero since u c = u T + u p . As the result, 
the mean velocity profile U fu c is significantly affected by the adverse pressure gradient (see the 
left of Figure 6). For flows with favorable pressure gradients, as we have discussed before that 
the ratio u p /u c for boundary layer flows could be larger than the value for channel or pipe flows, 
0.044, hence, the effect of favorable pressure gradients could be also significant. This is shown in 
the right of Figure 6. 




Figure 6: Effect of pressure gradients on the surface flow 


6 Conclusion 

A generalized wall function has been derived, which accounts for the effect of various adverse and 
favorable pressure gradients and is valid from the solid wall up all the way to the inertial sublayer. 
We have demonstrated that the effect of pressure gradients on surface flows are very significant, 
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especially for the flows with adverse pressure gradients. The traditional wall function must be 
replaced by the new generalized wall function for flows with considerable pressure gradients. We 
expect that this generalized wall function will be useful for analytical studies of turbulent surface 
flow properties as well as for CFD applications in providing appropriate boundary conditions for 
high Reynolds number turbulent complex flows. 
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